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SUMMARY 

Tracking control of the Space Station Freedom solar array beta gimbals is investigated. Of particu- 
lar interest is the issue of control in the presence of uncertainty in gimbal friction parameters. Sensitivity 
functions have been incorporated into the feedback loop to desensitize the gimbal control law to param- 
eter variations. Simulation results indicated that one such sensitivity function improves the closed-loop 
performance of the gimbals in the presence of unexpected friction parameter dispersions. 


INTRODUCTION 

The problem of control of mechanical joints in the presence of friction has been of considerable 
interest for years. Most large mechanical systems require actuation, and friction is inevitably present in 
such actuators. A principal issue in pointing control is the modeling of friction in gimbaling joints (ref. 5) 
Rockwell International’s Rocketdyne Division, with responsibility for the Space Station Freedom’s solar 
array beta gimbals, has investigated the dynamics of the beta gimbals (refs. 2 and 6) in the presence of 
four types of friction: Coulomb, Dahl, static and viscous. 

Both Coulomb and Dahl friction are nonlinear discontinuous functions of gimbal velocity, while 
static friction represents a dead zone. In this study, a quadratic cost criterion is applied to control design 
for the beta gimbals. The cost criterion penalizes tracking deviations and energy input to the motor. 

A term may be included in such a cost criterion which takes into account uncertainty in friction 
parameters, with the use of a sensitivity model in the feedback loop. The sensitivity model derives the . 
partial derivative of the state with respect to uncertain parameters, and provides a measure of change in 
state with parameter dispersions. 

Friction models incorporated in the gimbal dynamics are based on theoretical and experimental 
analysis performed by NASA Lewis Research Center and by Rocketdyne. Viscous friction is linear, and 
therefore simple to model. Coulomb friction is a constant torque that opposes the direction of gimbal 
velocity. Dahl friction (ref. 3) has a first-order dependence on the direction of gimbal velocity and is also 
highly nonlinear. Static friction is considered to be an algebraic inequality, not present directly in the 
dynamics. In the case of each form of friction, the tools of sensitivity analysis are applicable. 


The procedure followed here is as follows: (l) Using the known beta gimbal dynamics, sensitivity 
models are derived for each uncertain parameter. (2) A linear control law is applied which incorporates a 
penalty on the sensitivity functions. (3) Simulation results are used to assess and reduce the sensitivity 
models. 

NOMENCLATURE 


Name Description Units 


/3 gimbal angle radians 

I M motor current amps 

J A array axis moment of inertia in-Ibf-sec 

L m motor inductance H 

R m motor resistance 0 

KgEMF motor back-EMF constant volt-sec 

K t motor current-to-torque constant in-lbf/ amps 

Ky viscous friction in-lbf-sec 

T a applied gimbal torque in-lbf 

Tq Coulomb friction in-lbf 

T D Dahl friction in-lbf 

a slope of Dahl friction model in-lbf 

T dl Dahl friction limit in-lbf 

T s static friction in-lbf 

Tp total friction torque in-lbf 

T M motor torque in-lbf 

Ty viscous friction in-lbf 

V M motor voltage volts 


SENSITIVITY MODELS 

The known beta gimbal dynamics are shown in block diagram form in figure 1. The block diagram 
(refs. 2 and 6) may be used to prepare the model’s governing equations: 
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The Dahl friction term (ref. 3) represents bearing friction and is hysteretic. It can be represented 
more simply than in equation (l), although the above representation is widely used in ball bearing 
simulation studies. 

Equation (1) can be represented by a state-space description. In the case when the gimbals are moving, 

X = flX,a) + g(X)u 

where X = [ l u P T D ] T ^ 2 ) 

« - Vm 

« - [* n T c T DL f 

and f(X,a) and g(X) are obtained directly from equation (1), the “else condition. 

The sensitivity functions are defined as the partial derivative of the state vector with respect to the 
parameter of interest (ref. 4) . They are generated as follows: 


3 


dX 

da, 


M 

dX 




T 2 (X,«) 


M. M. + -V. (,* = 1,23) 

ax da, da, 


BEMF 

Ly 

1 


A 

0 


0 
0 

a, 1 

o — - 

h J A 

0 Tj(X,a) x 2 (X,«) 


1 - — sgn(xj I sgn[l - — 5 ^ti(Xj)] 

B 3 J ®3 

x A sgn(xj x A , 

-2a [1 a^Kxj)] — sgn[l x. 


JL = 

da, 

J L = 

da 2 

M. 

da 3 


[0 0 -X 3 Of 

[ 0 0 -sgn(x j) 0 ] r 


0 0 0 -t 2 (X,*)— 


( 3 ) 


This is linear in the sensitivity functions although nonlinear in X. The first of equation (3) is a linear, 
though time-varying, differential equation. 

The sensitivity functions are generally of interest in a neighborhood of some nominal parameter 
values. This implies that for all instances of a j the sensitivity model would be known. This is useful for 
small deviations in the friction parameters from their assumed nominal values. “Small deviations” Me 
considered here to be within 25 percent of nominal. For larger dispersions, second-order sensitivity 
functions may be needed (ref. 4), and those are a straightforward extension of equation (3). 

Having obtained the sensitivity model, to provide a linear control law that penalizes sensitivity merely 
involves a feedback gain matrix on the sensitivity, as in figure 2. Simulations of the model shown in figure 2 
were performed, in order to assess gimbal control performance improvement due to sensitivity feedback. 


SIMULATION RESULTS 

A truth model was prepared in FORTRAN, containing equation (2). A PID control law was designed 
to enable the gimbals to perform well based on nominal values of friction parameters, chosen as follows: 
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Value 

K vs 

81.4 

T c 

6.0 

t dl 

24.0 

T s 

12.0 

<7 

900.0 


The PID control law had the following gains: 


Gain 

Value 

Units 

proportional 

7430 

volts/rad 

integral 


volts/rad-sec 

derivative 

17530 

volts-sec/rad 


To incorporate all sensitivity functions would be equivalent to implementing a 12th-order filter in 
the feedback path, requiring more constants than are available to the gimbal controllers (ref. 7) in flight. 
The analysis was therefore used to select only those sensitivity functions which offered the greatest bene- 
fit to dynamical response. 

Gimbal position sensitivity functions offered little impact. Dahl friction and motor current states 
were found to be less important than gimbal velocity. Further, Coulomb and Dahl friction sensitivity 
parameters require differentiation across discontinuities, which can be implemented but should be avoided 
in this case due to resource limitations. For these reasons, only one sensitivity function was implemented 
in the control law: that of gimbal velocity (x 3 ) with respect to viscous friction (oj). 

Simulation results show that feedback of this sensitivity function will not only reduce control sensi- 
tivity to variations in viscous friction, but Coulomb friction variations as well. The function will not 
improve performance under Dahl or static friction dispersions, but will not degrade performance either. 

For a gain of 50 V-in-lbf-sec 2 /rad on this function, the following is a comparison of performance (in 
terms of percent overshoot and settling time) for a unit step command to gimbal angle, for cases with 
(W) and without (W/O) the sensitivity function feedback, for a 20 percent dispersion in K v : 


Indicator 

Nominal 

20 percent viscous dispersion 

W 

W/O 

increase 

decrease | 

w 

W/O 

w 

W/O 

overshoot 







(percent) 

10 

10 

8 

8 

9 

9 

settling 







time (s) 

28 

35 

28 

35 

33 

42 


The following is a comparison of performance for a unit step command to gimbal angle, for cases 
with and without the sensitivity function feedback, for a 20 percent dispersion in T^,. 
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Indicator 

Nominal 

20 percent Coulomb dispersion 

w 

W/O 

increase 

decrease jj 

w 

W/O 

w 

W/O 

overshoot 







(percent) 

10 

10 

8 

8 

8 

8 

settling 







time (s) 

28 

35 

30 

30 

30 

90 


As can be seen from these two performance summaries, the only case in which no measurable performance 
improvement can be seen is that of a 20 percent increase in Tq. Even in that case, there is no perform- 
ance degradation. 

In each case, the sensitivity feedback gain is maintained at 50. This gain can be changed as the 
designers gain confidence in friction parameter values. If, for instance, there is a tolerance of 5 percent or 
less expected in Coulomb and viscous friction parameters, then this feedback term will not improve 
performance much. If, however, the tolerance is as high as 25 percent, then the introduction of this term 
to the gimbal control law will improve performance a good deal. 

CONCLUSIONS 

The addition of a single sensitivity function to the feedback control law, along with an existing 
PID controller, improves settling time without increasing overshoot in step response. Beta gimbal dynam- 
ics were simulated both with and without this sensitivity function. Simulation results showed improved 
performance for parameter variations in viscous, Coulomb and Dahl friction from nominal values. 

The sensitivity function, that of gimbal velocity with respect to viscous friction coefficient, is a 
recommended addition to the beta gimbal control law algorithm for Space Station Freedom solar array 
pointing control. 

Further work, if any, should consider the inclusion of other degrees of freedom, such as station 
body rotations and array flexible modes. 
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Figure 1 .—Block dia&am of beta gimbai control and dynamics without sensitivity. 



Figure 2.— Block diagram with sensitivity feedback of gimbai velocity. 
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